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Abstract. We consider eigenvalues of a product of n non-Hermitian, in- 
dependent random matrices. Each matrix in this product is of size N x N 
with independent standard complex Gaussian variables. The eigenvalues of 
such a product form a determinantal point process on the complex plane 
(Akemann and Burda [T|), which can be understood as a generalization of 
the finite Ginibre ensemble. As N — > oo, a generalized infinite Ginibre 
ensemble arises. We show that the set of absolute values of the points of 
this determinantal process has the same distribution as {-Rj" , R^ > • • •}> 

where are independent, and (^R^^j is distributed as the product of 

n independent Gamma variables Gamma(fc, 1). This enables us to find the 
asymptotics for the hole probabilities, i.e. for the probabilities of the events 
that there are no points of the process in a disc of radius r with its center at 
0, as r — ¥ oo. In addition, we solve the relevant overcrowding problem: we 
derive an asymptotic formula for the probability that there are more than 
m points of the process in a fixed disk of radius r with its center at 0, as 
m — >■ oo. 



1. Introduction 

Products of random matrices are used in different areas of research. For ex- 
ample, the book by Crisanti, Paladin and Vulpiani [9] describes applications of 
products of random matrices in statistical mechanics of disordered systems, lo- 
calization, wave propagation in random media, and chaotic dynamical systems. 
For an application of such products to the study of compositions of random 
quantum operations we refer the reader to the paper by Roga, Smaczyhski 
and Zyczkowski [30J. A paper by Osborn [27] considers products of random 
matrices in the context of Quantum Chromodynamics. For applications to 
Schrodinger operators see the book by Bougerol and Lacroix [5]. 
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In the 1960's and 70's different fundamental probabilistic results on prod- 
ucts of random matrices were obtained. In particular, the asymptotic behavior 
of products of independent random matrices was investigated in the work by 
Furstenberg and Kesten [20]. However, in the 1960's and 70's spectral aspects 
of products of random matrices did not attract any serious attention of math- 
ematicians working in the field. Only recently a number of works appeared in 
which eigenvalue distributions of products of random matrices, in the limit of 
large matrices, were considered (see, for example, Burda, Nowak, Jarosz, Livan 
and Swiech [HIE], Burda, Janik, and Waclaw [8], Gotze and Tikhomirov [H] . 
Penson and Zyczkowski [29J, O'Rourke and Soshnikov [32J, Forrester [IH]), and 
where products of random matrices were studied by usual methods of Random 
Matrix Theory. We refer the reader to the books by Anderson, Guionnet and 
Zeitouni [3], Deift [10], Forrester [17], and Pastur and Shcherbina [28] for an 
introduction to Random Matrix Theory, and for the description of its basic 
methods and results. 

In this article we concentrate on radial distributions of eigenvalues of prod- 
ucts of complex non-Hermitian independent random matrices. For a properly 
normalized product of complex non-Hermitian independent random matrices 
O'Rourke and Soshnikov [32] showed (under certain assumptions on the en- 
tries of the random matrices) that the empirical spectral distribution of the 
eigenvalues converges to the limiting distribution, which is a power of the cir- 
cular law. Forrester [19] derived a formula for the Lyapunov exponents for a 
product of complex Gaussian matrices. The starting point of the present re- 
search is the result obtained in Akemann and Burda pQ. They considered the 
product of n complex non-Hermitian, independent random matrices, each of 
size N x N with independent identically distributed Gaussian entries (Ginibre 
matrices). It was shown that the eigenvalues of such a product form a complex 
determinantal point process which can be understood as a generalization of the 
classical Ginibre ensemble. It is the aim of the present paper to study in detail 
this complex determinantal process (which in this paper is called the gener- 
alized finite- N Ginibre ensemble with parameter n), and its infinite analogue 
(which in this paper is called hereafter the generalized infinite Ginibre ensemble 
with parameter n). We show that the set of absolute values of the points of 
the generalized Ginibre ensemble has the same distribution as {b!{ \ • • •}, 

where are independent, and (jrfj^j is distributed as the product of n 

independent Gamma variables Gamma(&,l), see Theorem 13.11 and Theorem 
13.21 This enables us to find the asymptotics for the hole probabilities, i.e. for 
the probabilities of the events that there are no points of the process in a disc 
of radius r with its center at 0, as r — > oo, both for the generalized finite- N 
and for the generalized infinite Ginibre ensembles (Theorem 13. 5j) . In addition, 
we solve the relevant overcrowding problem: we derive an asymptotic formula 
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for the probability of the event that there are more than m points of the gener- 
alized infinite Ginibre ensemble in a fixed disk of radius r with its center at 0, 
as m — > oo, see Theorem 13.61 In proving these Theorems we apply a technique 
similar to that developed in Kostlan [23], Krishnapur [21] and Hough, Krish- 
napur, Peres and Virag [2TJ for the case of the classical Ginibre ensemble, and 
for the case of random analytic functions. 

This paper is organized as follows. In Section [2] we describe the result ob- 
tained in Akemann and Burda [TJ, and define explicitly the relevant determi- 
nantal process. In Section 3 we state main results of this paper. Theorem 
13.11 and Theorem 13.21 describe the distribution of absolute values of the points 
of the generalized Ginibre ensembles, Theorem 13.41 gives an exact formula for 
the hole probabilities, Theorem 13.51 concerns the rate of the decay of the hole 
probabilities, and Theorem 13.61 solves the relevant overcrowding problem. The 
rest of the paper is devoted to proofs of these results. 

Acknowledgements. Part of this research was conducted during ZIF research 
program " Stochastic Dynamics: Mathematical Theory and Applications" . It is 
our pleasure to thank the Center for Interdisciplinary Research (ZIF) of Bielfeld 
University for hospitality, and the organizers of the ZIF Research Group 2012 
"Stochastic Dynamics: Mathematical Theory and Applications" for the stim- 
ulating and encouraging environment they created at the program. 

2. Products of random matrices and generalized Ginibre 

ensembles 

In this article we consider the product P n of n independent random matrices, 



Each matrix Xj, j — 1, . . . , n, is of size N x N, and with i.i.d standard complex 
Gaussian entries. Let zi,. . .,zn be the eigenvalues of P n . We study some statis- 
tical properties of the distribution of the eigenvalues z\ y . .,Zn in the complex 
plane, both for finite and for large (N — > oo) matrices. More explicitly, the 
starting point of the present work is the following result obtained recently by 
Akemann and Burda [TJ. 

Assume that \zi\ < ... < \z N \. Then the joint density of (2^)4=1,. ..,tv with 
respect to Lebesgue measure on is given by 



n 



X\ X2 ... A, 



) 



n 



N 



(2.1) p%\ Zl ,...,z N ) 



1 
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k=l l<i<j<N 
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Here G : 



n,0 
0,n 



0, 0, 



I stands for Meijer's G-function with suitable 



choice of parameters. For Meijer's G-functions we adopt the same notation and 
definitions as in Luke [2E], Gradshtein and Ryzhik [T5]. Namely, the Meijer 

ai ' a2 ' ' ' ' ' ap lk <],.(m,->(i 
77 1 I lb LlcllIlcLL do 

h, 02, • Oa 



G-function G™' q n I x 



fim,n I „ 



bi, b 2 , 



1 

2vn 



UT=inb,-s)UU T (l-a J + s) 



c 



it 



j=m+ 



Here an empty product is interpreted as unity, < m < q, < n < p, and 
the parameters {a k } (k = 1, . . . ,p) and {bj} (j = 1, . . . , m) are such that no 
pole of Y(bj — s) coincides with any pole of T(l — a k + s). We assume that 
z £ C \ {0}. The contour of integration G goes from —zoo to -Moo so that 
all poles T(bj — s), j = 1, . . . , m, lie to the right of the path, and all poles of 
r(l — a k + s), k = 1, . . . , n, lie to the left of the path. If p = 0, then n = 0, and 

we write the corresponding Meijer G-function as bi, b 2 , b, 

In particular, we have 

c+ioo 

Gjftf |0, 0, . . . , 0) = JL J r s T n (s)ds (t > 0, c> 0). 

c—ioo 

This integral can be evaluated (see Springer and Thompson [33] , Lomnicki [25] ) 
by contour integration in the form of an infinite series: 

oo 

Gg(t|0,0,...,0) = ^#(t,n,j), 

3=0 

where R(t, n, j) is the residue of the integrand at the nth-order pole 

s = -j (.7=0,1,...), 

i.e. 



i Jn-l 



n — 



Specifically, we find 



G^(i|0,0,...,0) = ^(-1 



tP 



3=0 



and 



V 



G 2 '° 2 (t\0, 0, . . . , 0) = ^ — ^(- log* + 2V(j + 1)) = 2K (2Vi). 



i=o 
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Here is the Euler psi function and Kq(.) is the modified Bessel function of 
the second kind of zero order, see, for example, Erdelyi [T3] . 

Equation (12. ip implies that the eigenvalues of P n form a determinantal point 
process^ on the complex plane with kernel 



N-l ( f\k 
fc=0 v ' ; 

with respect to the background measure ^w n (z)dm(z). Here dm(z) denotes 
the Lebesgue measure on the complex plane. The fact that ■^w n (z)dm(z) is 
indeed a probability measure for any positive integer n can be checked using 
the formula 



oo 



lu— 1 r-im,n 
1 ^p.g 



di, d2, . . • , a p 



q 



dt 



(2.2) 



]1 T(b s + u)U r(l - a, - u) 
j=i i^i 

ft Y{l-b 3 -v) f[ V{a J + v) 

j=m+l j=n+l 



see Luke [26], Section 5.6.7. 

We will refer to the determinantal point process formed by the eigenvalues 
of P n as to the generalized finite-N Ginibre ensemble with parameter n. The 
reason is that once there is only one matrix in the product (n = 1), we have 

w 1 (z) = G 1 %z\ 2 \0,0,...,0) = e-^, 
and equation (12. ip turns into 

1 N 



^ 11^=1 r (^) fc=l l<i<i<Ar 



Therefore the determinantal point process formed by eigenvalues of P n reduces 
to the classical Ginibre ensemble at finite N, i.e. to the determinant point 
process with the kernel 



*r i} (*,o=E 

k=0 



k\ 



with respect to the background measure -e - '*' dm(z). 

If iV = oo, then we call the corresponding determinantal point process the 
generalized infinite Ginibre ensemble with parameter n. 



1 For a background on determinantal point processes we refer the reader to survey articles 
by Borodin 0], and by Hough, Krishnapur, Peres, and Virag [52]. 
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3. Statement of results 

3.1. The distribution of the moduli of eigenvalues. Our first result con- 
cerns the distribution of the moduli of the eigenvalues of P n . Recall that gamma 
variables Gamma(A;, 1) are those having the following density function 

(3.1) «?><*> = { o**"^' HI 

Theorem 3.1. Let P n = X\Xi . . .X n be a product of n independent random 
matrices. Each matrix Xj, j = l,...,n, is of size N x N, and with i.i.d. 
standard complex Gaussian entries. The set of absolute values of eigenvalues 

of P n has the same distribution as the set ^R^\ . . . , R$\ , where the 
random variables R[ , . . . , R^ are independent, and for each k, 1 < k < 
N, the random variable (r^} has the same distribution as the product of n 
independent and identically distributed gamma variables Gamma(fc, 1). 

The next Theorem gives even more explicit information on the set of random 
variables { R^ , R^ , . . . , R ( f 



Theorem 3.2. The random variable (r^^\ has the density function given 
by the formula 



k-1, k-1, k- 1 ) , x>0, 



(3.2) gf\x) = { (r(fc))"^o,n 

0, x < 0. 

Remarks. 

(a) Theorem 13.11 generalizes the result obtained by Kostlan [23] in the case of 
the classical Ginibre ensemble (n = 1). 

(b) The generalized finite- N Ginibre ensemble with parameter n is a determi- 
nantal process on C with the kernel of the form K(z,£) = ^2,Ck{z^ ) k (where 

k 

Ck are some coefficients) with respect to a radially symmetric measure. It is 
a known general fact (see Hough, Krishnapur, Peres and Virag |2T], Section 
4.7) that the set of absolute values of the points for such processes has the 
same distribution as a set of independent random variables. Theorem 13.11 and 
Theorem 13.21 describe this set of independent random variables explicitly. 

(c) Theorem 13.21 is closely related to the following result on the distribution of 
a product of independent gamma variables (see Springer and Thompson 
Section 3). 
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Proposition 3.3. Let X\, x%, . . ., x n be n independent gamma variables having 
density functions 

1 x^e-**, x k >0, 



10, x k < 0, 

where b k > 0, 1 < k < n. Then the probability density function g(z) of 
the product z = x\Xi . . . x n is Meijer's G -function multiplied by a normalizing 
constant, i.e. 

g(z) = Go'°(z\h - 1, b 2 - 1, . . . , b n - 1). 

n r(6,) 

i=l 

3.2. An exact formula for the hole probabilities. Denote by Nqq{t\ N) 
the number of points of the generalized finite- N Ginibre ensemble with param- 
eter n in the disk of radius r with its center at 0. Alternatively, Nqq(t\ N) can 
be understood as the number of eigenvalues of the random matrix P n in the disk 

of radius r with its center at 0. By the hole probability Prob jA/^(r; N) — oj 

we mean the probability of an event that there are no points of the generalized 
finite-iV Ginibre ensemble with parameter n in the disk of radius r with its 
center at 0. 

Theorem 3.4. The hole probability Prob jA/^^r; N) = j for the generalized 
finite-N Ginibre ensemble with parameter n can be written as 



Prob{A/S(r;A0=0}=;Q 

k=l 

Remarks. 



n n+l,0 I 9 

N G l n+1 [ r 



1 

, J\i j • * * J 



(T(k)T 



(a) Note that there is a convenient integral representation for the Meijer G- 
function in the formula above, namely 

c+ioo 



(b) Since 



1 

0, k, 



)=± f r ^(T(k + s)r^, c>0. 
J 2m J s 



1 

0, k 



T(k,r 2 ) = / e'H^dt 



we see that once n = 1 the formula in the statement of Theorem 13.41 reduces 
to 

Prob^r\r;N) = 0} = f[ T -^l. 

k=i ^ ' 



cS 
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The formula just written above for Prob |a/^ ^(r; N) = j is well known, 
see Grobe, Haake, and Sommers [16 1, Forrester 1181. 



3.3. The decay of the hole probabilities. A basic quantity of interest is 
the decay of the hole probability as r — > oo. We investigate the decay of 
the hole probabilities both for the generalized finite- N Ginibre ensemble with 
parameter n (the case of products of n random matrices, each of which is of 
size N), and for the infinite generalized Ginibre ensemble (the case of product 
of n infinite random matrices). 

Theorem 3.5. (A) (The case of products of random matrices of finite size 
N). 

As r oo, 



Prob 



{A/S(r;iV)=0} 



(ra— 1)N 

exp | -niVr» + N In j log(r) I (l + O (V 



M N 



ni n( r (fc)) n 

k=l 



This implies that for the product of n matrices of finite size N we have 

1 



lim ( — lo£ 



Prob { (r;N) = 



-nN. 



(B) (The case of products of infinite random matrices). 
The following limiting relation holds true 

lim ( 4- log [Prob [U^Xir] N = oo)\ = ol ] 

r^oo \ r - LI J \ J 



n 
4' 



Remarks. 

(a) For the classical (infinite) Ginibre ensemble we have n — 1, and Theorem 
13.51 (B) sa ys that 

1 Prob {logA4 n G =1) (r; N = oo) = o} -+ -1 

as r — > oo. This asymptotic result for the classical Ginibre ensemble is well 
known, see Grobe, Haake, and Sommers [16J, Forrester [18J, Hough, Krish- 
napur, Peres, and Virag [21], Akemann, Phillips, and Shifrin [2] for different 
proofs and related results. 

(b) The result of Theorem 13.51 (B) can be compared with the decay of hole 
probabilities for the zeros of the Gaussian analytic function, 
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where a n are i.i.d. standard complex Gaussian random variables. Namely, it 
was proved by Sodin and Tsirelson [21] that the hole probability for the zeros 
decays like exp{— cr 4 }. Theorem 13.51 says that for the generalized infinite Gini- 
bre ensemble with parameter n the hole probability decays like exp{— Cr~}. 
Thus we have the same behavior only for n = 1. 

3.4. Overcrowding. Consider a disk with a fixed radius r > 0, and with its 
center at 0. Recall that Nqq(t\ N = oo) denotes the number of points of the 
generalized infinite Ginibre ensemble with parameter n in this disk. The prob- 
lem is to estimate the probability of the event that in this disc there are more 

than m points of the ensemble, i.e. to estimate Prob |A/"(^](r; N = oo) > m j. 

We are especially interested in the decay of this probability, Prob vA/^Cr; iV = oo) > m 
as m — ¥ oo. 

In this article we prove the following Theorem. 

Theorem 3.6. Let M^q{t\ N = oo) be the number of points of the infinite 
generalized Ginibre ensemble in the disk of radius r around 0. Then for a fixed 
r > 

Prob |A/"(^](r; N = oo) > m j = exp | — -nm 2 log(m)(l + o(l)) j , 

as m —¥ oo. 
Remarks. 

(a) Theorem 13.61 is a generalization of the result obtained by Krishnapur (see 
Krishnapur [21], Section 2.1) for the classical infinite Ginibre ensemble (the 
case corresponding to n — 1). 

(b) In the context of zeros of Gaussian analytic function the overcrowding prob- 
lem was formulated by Yuval Peres, and was studied in detail by Krishnapur 
in [21]. It was shown that the probability of the event that in the disk with a 
fixed radius r around there are more than m zeros of the Gaussian analytic 

function decays in the same way as Prob Y^GG~ 1 \ r: > N = oo) > m j. 

4. Proofs of Theorem 13.11 and 13.21 

Let ri, . . . ,rjv be the moduli of the eigenvalues z±, . . . ,zn of P n , i.e. r\ = \zx\, 
r 2 = l^l, • • •, fN = \zn\- Thus ri < . . . < r^v. We want to find the joint density 
of {ri) i=l! ... ;N . 

Proposition 4.1. The joint density of (r^i^.^ is given by 

(27T) 



Per [r^ ^^H^ifa)- 
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Proof. Let z±, . . . , be the eigenvalues of P n . The joint density of (2^)4=1,...,^ 
is given by formula (12. ip . Set 



We have 



n 

l<i<j<N 



Zi = ri e M \ i = l,..., N. 



A 



Z'i Zj I 



e (-i) sgn(ff) n 

<reS(AT) i=l 



E 

a,a'eS(N) 



N 



A 



'_ 1 }sga(<r)+egn(«T / ) JJ r jC?')-l e i(o-(j)-l)% JJ r . CT '( ft )- 1 p -i(^'(fe)-l)^ 



r fc e 



fe=i 



This gives 



2?r 2tt 



iV 



n 

l<i<j<N 



Zi - Z 



: j \ 2 d6 1 ...d6 N = {2n) N E II 



,2a(j)-2 



o-eS(N) j=i 

Therefore the joint density of (r i )j = i v .. i Ar is 

N N 

U w n(r k ) e n 



(2tt 



1 



f2 7 r) iV 



2<t(j)-2 



fc=l 



ffeS^AT) j=l 



N 

n 



A" 



T"n£*r(fc))» 



p er K 2j 1 ]£=iII u,n ( r ^ 



□ 



Proposition 14.11 implies that the random variables yi = rj, 1/2 = r§, . . ., 
Un = r 2 N have the joint density given by 



per 



-I AT 



_7r»~ 1 (r0')) ; 

Observe that formula (12. 2p implies 



J i,j=l 



y*- 1 w n {y/y)dy = ir- 1 r n (j). 



Therefore the vector of squares of absolute values of eigenvalues of P n (in 
uniform order) has the density 



M per 



ef\m) 



1 N 
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where the functions Qj(y), 1 < j < N, are probability density functions 
defined by 

3 \ 0, 2/ < 0. 

Since 

y^wniVv) = iT^G/lo, o, . . . , o) = G${ y \j - 1, j - 1, . . . , j - 1), 

we conclude that gf^iv) can be rewritten as in equation (13.21) . 

To complete the proof of Theorem 13.21 we use the following well known fact 
(see, for example, Kostlan [23], Lemma 1.5). Assume we are given an iV-tuplet 
of independent random variables (A), 1 < i < N, with the corresponding 
densities (ft), 1 < i < N. Define a new iV-tuplet of random variables, (Bi), 
1 < i < N, as a random permutation of the vector (A), 1 < i < N (these 
random permutations are equal to each other in probability). Then the joint 
density of the random vector (Bi), 1 < i < N, is 

^per re- 
considering squares of moduli of unordered eigenvalues of P n as random vari- 
ables (Bi), 1 < i < N, we obtain the statement of Theorem 13.21 

Theorem 13.11 follows from Theorem 13.21 and from the result by Springer 
and Thompson [34J on the distribution of a product of independent gamma 
variables, see Proposition 13.31 Namely, Proposition 13.31 and Theorem 13 . 2 1 imply 

I n\ 2 



that each random variable (R\ j , 1 < i < N, has the the same distribution 

as the product of n identically distributed and independent gamma variables 
having density function (I3.ip . Theorem 13.11 is proved. □ 



5. Proof of Theorem 13.41 
From Theorem 13.11 we conclude that 

N , 

Prob {jV§g(r; N) = o} = Prob [R k 

k=l ^ 

where the random variables Rx, R% , ■ ■ ■, R^ are those introduced in the 
statement of Theorem 13.11 By Theorem 13.21 



1 

(my 



Prob {(R^ ) >r 2 \ = 7F77X ^ / G n ° n (x\k -l,...,k- l)dx. 
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The last integral can be computed explicitly using the formula 



a ll • ■ ■ ! 

p-cj, 6l, 



see Gradshtein and Ryzhik [15J, 7.811.3. This gives 



Prob \ (Rf) 2 > r 





-1, jfe-1, fc-1 



Since 



ai, 
61 , 



"p,5 1 ^ 



(r(*))» 



ai + cr, . . . , a p + cr 
61 + a, . . . , 6q + a 



(see, for example, Luke [26], Section 5.4) we can rewrite the expression for 



ProM >r 2 > as 



ProM (4 n) T >r 2 



; k ; . . . . 



(r(fc)) n 



and the result of Theorem 13.41 follows. 



□ 



6. Proof of Theorem 13.51 

6.1. Proof of the asymptotic formula for the hole probability for the 
generalized finite- N Ginibre ensemble. We use Theorem 13.4} which ex- 
presses the hole probability Prob jA/^gfV; N) = Oj in terms of the Meijer G- 

functions. The following asymptotic formula holds true (see Luke j2S] , Section 
5.7) 



G&i, (I2, 

h, b 2 , 



, dp 

, bg 



(2^)(-D/2 



a 



1/2 



exp 



k z 



-k/a 



k=0 



where 

\z\ — > 00, I argz| < (a + e)n — 5, 5 > 0. 
In the formula above the parameters a and e are defined by 



and 



a = q-p, 



e = - if (7=1, e = 1 if a > 1. 
2 
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The parameter 9 is defined by the formula 



f 1 Q P } 

I k=i fc=i J 



Finally, the coefficients M^'s are independent of z and can be found explicitly. 
In particular, M = 1. 

In our case p = 1, q = n + 1, a = n, and z = r 2 . It is not hard to find the 
parameter 9 as well. The result is 



2 2n 



This gives 



G n+ifi 2 
^l.n+l 1 ' 



1 \ (2ti)— r 2^1 2fc i_i .( i 

— j — exp < — nr « > r 



, A? , . . . , fv 



«2 



1 + (r~- 



as r — >■ oo. We insert this asymptotic expression into the formula in the 
statement of Theorem 13.41 The statement of Theorem 13.51 (A) follows imme- 
diately. □ 

6.2. Proof of the asymptotic formula for the hole probability for the 
generalized infinite Ginibre ensemble. 

6.2.1. An upper bound for the hole probability. To estimate the hole probabil- 
ities we use the following standard fact (called the Markov inequality). 

Proposition 6.1. Suppose : K — > R is a positive valued function, and let A 
be a Borel subset o/R. Then 

inf {<p(y) :yeA}- Prob {X G A} < E(p(X). 

Proof. See, for example, Durrett [12], Section 1.6, Theorem 1.6.4. □ 

Proposition 6.2. We have 

Prob |A^(r; iV = oo) = o} < exp [~r$ + O (log(r))} , 
as r — > oo. 

Proof. Let a > 0, A = (r 2 , oo), and ip(x) = x a . Then by the Markov inequality 
(Proposition 16.11) we have 

oo 

ft<*cr >° n (t\k-i,k-i,...,k-i)dt 

Prob<! IR^) >r 2 \ < 



* ) ( - (r(fc))"(r 2 ^ 
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(we have used the fact that the random variable yR^j has the density func- 
tion given by formula ( 13.21) . see Theorem I3.2[) . By formula (12. 2\\ 

oo 

J t a G^(t\k - 1, k - 1, . . . , k - l)dt = (T(k + a)) n . 



Therefore, 



Prob^(^ 2 >r4< (r(fc + Q))? 



k I I - (r(A;))«(r 2 ^ 



Next we use the following well known inequality (see, for example, Digital 
Library of Mathematical Functions [TTj, § 5.6) 
(6.1) ' 

(27r)^exp \-z+ (z ) log si < T(z) < (27r)^exp \-z + (z - -)\og{z) + ! 



as z > 1, to obtain 



Prob<! [R^Y >r 2 

< exp < -na - alog(r 2 ) +n ( k — \ + a J log f 1 + ^) + anlog(k) + W 

[ \ 2 J \ ks 12{k + a) 

2 

Set a = r» — k. Then we can rewrite the inequality above as 

f/ ( n )\ 2 2 1 f ,2 . n k k n 

Prob < I R y k J > r > < exp < —n(r « — fe) H — log — 5 — nfc log — + 



2 7~n }"n 12r 



This gives 



2_ 

JI Prob (R^) > r 2 \ < \{ Prob (r { ^ > r 2 \ 

k=l ^ ' k=l ^ ' 

f 2 2. 
-nr-n + n + - ^ log ( I - n ^ Hog ( \ + - 
fc=i ' fc=i 



where r ™ is considered as an integer (this assumption should not affect our esti- 
mate). The sums in the exponent can be estimated using the Euler-MacLaurin 
sum formula. We write it in the form 

l L 

(6-2) fl W = f f(t)dt + \ (f(L) + /(l)) + O U'{L)) . 

fc=i { 
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This formula gives 



2 
r 'J 



E lo f 

k=l 



k 

2 

T n 



and 



2 
r n 



2 1 / 2 

-r n + 1-| — log(r) + O I r~« 
n V 



as r — > oo. Using these estimates we find 

Prob {A/S(r; AT = oo) = o}=JJ Prob I (r^A > r 2 I 

k=i ^ ' 

< exp | _ ^ r ™ + O (log(r))| , 

as r — > oo. 

6.2.2. A lower bound for the hole probability. 
Proposition 6.3. We have 

Prob {-A/^r; N = oo) = oj > exp {-^" + 0(r% log(r)) } , 
as r — > oo. 

Proof. It is known that 



(6.3) 



n Pr ° b 



k=i 



R 



(n=l) 



> r 2 \ > exp < — r 4 + 0(r 2 log(r)) 
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□ 



see, for example, Hough, Krishnapur, Peres and Virag [2TJ, Section 7.2. Con- 
sider the set of the random variables {R[ , i?2 , • • •}• We know that the ran- 
dom variables are independent, and yR^j has the same distribution 
as the product of n independent and identically distributed gamma variables 
Gamma(fc, 1). In particular, the random variable (^R^~^ J is itself the gamma 

variable Gamma(/c, 1). We conclude that the random variable [R^ ) has the 



same distribution as the random variable 
implies 



R 



(n=l) 



. This immediately 



(6.4) 



Prob <^ R 



>(») 



> 



Prob <^ R 



>(n=l) 



> r™ 
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Using inequalities ( 16. 3ft and ( 16. 4 j) we find 

Prob [jV^r; W = oo) = 0} = JJ Prob j (4*°) > r 2 I 

> exp | — — r™ + C(r" log(rj) 

as r — )• oo. □ 

6.2.3. Proof of Theorem \ 3.5\ (B). The statement of Theorem 13.51 (B) follows 
from Proposition 16.21 and Proposition 16.31 □ 

7. Proof of Theorem 13.61 

Recall that the set of absolute values of points of the generalized Ginibre 
ensemble has the same distribution as the set {R[ , . . .}, where R% are 

independent, and f-Rj^l has the same distribution as the product of n inde- 
pendent and identically distributed gamma variables having density function 
Gamma(A;, 1), see Theorem 13.11 This implies 

where \ l<i<k,l<j<n are i.i.d. exponential random variables with 
mean 1. Therefore we can write 

Prob | («<">) 2 <r 2 \> Prob {{|" + . . . + {<" < ri } ■ . . . ■ Prob {?«"' + . . . + < 'A } 

>np r ob(^»< 4}..... nprob 



For an exponential random variable £ with mean 1 we have 
This gives 



Prob {£ < x} > -, < a; < 1. 



{(*?') 



2 



Prob <^ < > — 



^ \ / J - 12*: 



r n 



and we obtain 



m ( 2 ^ 

Prob {-/V§3(r; N = oo) > mj > JJ Prob i (#[ n) ) < r 2 I 

fe=i ^ J 



^2fc „m(m+l) 

^ 11^^ = on ^( m+ D ex p i -»E Ho sW r ■ 

k=i * 2 



k=l 
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The Euler-MacLaurin formula (equation (16.2ft ) gives 

(7.1) 5>log(fc) = mlm 2 M og(m) - ^- + (tog(m)), 
fc=i 

asm-} oo. We use this estimate to get a lower bound for Prob j Nqq{t\ N = oo) > 
namely 

(7.2) Prob |A/^(r; TV = oo) > mj > exp j-^nm 2 log(m)(l + o(l)) J , 

as m — )■ oo. To obtain an upper bound for Prob jA/g^r; iV = oo) > mj we 

use the Markov inequality (Proposition 16. ip with A = {0,r 2 }, (p(x) = x 
a > 0. This gives 



as 

; 



oo 



J t- a GoZ(t\k — 1, fc — 1, . . . , fc — l)d* 

(7.3) ^b{(^<^}< M «i — : — ^ 

= r 2 r (r(fc-a)) n 
1 ; (r(ib)) n ' 

where we have used formula (12.21) . By inequality (16. ip 



Prob<! (i?l n) ) 2 <r 2 

< exp |na + a log(r 2 ) + n(k a) log(l — — ) — an log(fc) + 



2 ' &v fc 12(ife-a) 
Choosing a = A; — |, we obtain 



{7A) [( 1\ nl 

< exp j f fc - - J (n + log(r 2 ) - n log(fc)) + - > . 

In addition, by simple probabilistic arguments 



P]-ob^jV^(r;JV = oo; 



(7.5) 



> mj < Prob j^I (jR^) 2 < ^ m | 



+ Prob j (^i n) ) < r2 \ 

k=m?+l K ' 



(Here I (.) stands for the characteristic function of a set). The second term on 
the right hand side of the inequality above can be estimated as follows. By 
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inequality ( 17.4ft 

Prob|^ n) ) 2 <r 2 J < exp{-nA;log(/c)(l + o(l))} 



as fc — ?• oo. Therefore, 



__ Prob j(i?i n) ) 2 < r 2 | < exp{-nm 2 log(m 2 )(l + o(l))} 



k=m 2 +l 



as m — >■ oo. Now let us estimate the first term on the right hand side of 
inequality (I7.5p . We have 



Prob | __ I ((jZ<»>) 2 < r 2 ) > -} < ( ^ ) ft Prob { (itf < 

/ ?n 2 \ _____ 

Using I < m 2m , inequality ( I7.4p . and equation ( 17. ip we obtain 



Prob|^l((i4 n) ) 2 <r 2 ^j >i 

{nm m 1 
2mlog(m) H h ) ( n + l°g( r2 ) — nlog(fc)) 



nm 2 log(m) 
exp < (1 + o(l)) 



as m — > oo. The statement of Theorem 13.61 follows from inequalities (17. 2p and 

(T~). □ 
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